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Abstract 

We study the implications of lattice QCD determinations of the S-wave nucleon-nucleon scat- 
tering lengths at unphysical light quark masses. It is found that with the help of nuclear effective 
field theory (NEFT), not only the quark mass dependence of the effective range parameters, but 
also the leading quark mass dependence of all the low energy deuteron matrix elements can 
be obtained. The quark mass dependence of deuteron charge radius, magnetic moment, polar- 
izability and the deuteron photodisintegration cross section are shown based on the NPLQCD 
lattice calculation of the scattering lengths at 354 MeV pion mass and the NEFT power counting 
scheme of Beane, Kaplan and Vuorinen. Further improvement can be obtained by performing 
the lattice calculation at smaller quark masses. Our result can be used to constrain the time 
variation of isoscalar combination of u and d quark mass m^, to help the anthropic principle 
study to find the mg range which allows the existence of life, and to provide a weak test of the 
multiverse conjecture. 
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I. INTRODUCTION 



A very interesting aspect of lattice QCD (LQCD) calculations is that one can study 
the quark mass dependence of physical observables which are otherwise hard to measure 
with experiments. This information could be used to constrain the time variation jl-3| 
of quark masses in the evolution of the universe j4-7|. It could also shed light on how 
finely tuned the quark masses should be 8l4l2(| such that light nuclei can be synthesized 



13 



14j | and make the 



through the usual pathway of Big Bang Nucleosynthesis (BBN) 
familiar carbon based life forms possible. 

Much has been learned from the u and d quark mass (we will work in the isosym- 
metric limit m„ = = nig) dependence of the meson and single baryon observables 
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16| through lattice QCD (LQCD) |l7|, chiral perturbation theory (ChPT) [18|-l21| 



and experimental data. In principle, lattice QCD can map out all the rrig dependence 
for these observables. However, most of the calculations are done with m^'s larger than 
their physical values, because it requires more computing resources to work with smaller 



Fortunately, ChPT, which is an effective field theory (EFT) of QCD, can be used 



to described the dependence once the unknown parameters in the theory are fixed by 
either experiments or lattice data. 

In the multi-baryon sector, much progress has been made in LQCD in two nucleon 



22h26| . nucleon- hyperon 27|], triton 
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29| and a-particle [29i] systems (see [SOj for a 



brief review). However, for two nucleon systems, so far only the S-wave scattering lengths 
have been computed with 354 MeV or heavier pion mass m,r- (Note that the physical 



pion mass m: 



phys 



138 MeV, and there is a one-to-one correspondence between m-Tr and 

1/2 

niq, e.g. oc m^' as — )■ 0. So the m,r and niq dependence can be converted to 
each other.) Even so, as will be demonstrated in this work, this information is enough 
to determine the leading m^r dependence of all the low energy matrix elements involving 
deuterons. 

We will focus on processes with the typical momentum p <^ m-Tr, such that the pions 
can be taken as heavy particles and integrated out of the theory. This theory is known 



as pionless theory 
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L|-|34||. The information of the pion dynamics in the pionful theory 
is now encoded in the dependent couplings of the pionless theory. It is found that, 
all the leading dependence in deuteron matrix elements in the pionless theory can 
be computed using the pionful theory together with the dependence of the S-wave 



scattering lengths obtained from LQCD. Thus, once they are fixed at m^'*^*, their values 
at other pion masses are also known. 

Of course, one can still work with the pionful theory. The matching is a convenient 
but not necessary step to take. One advantage of working with the theory without pions 
is that once the dependence of the couplings are worked out, one can just perform the 
calculation in the pionless theory instead of the more complicated pionful theory. As an 
explicit example, we match the pionful theory based on Beane, Kaplan and Vuorinen's 



(BKV) 351] power counting scheme to a pionless theory. This allows the matching been 



done analytically. However, the method can be applied to other power schemes as well. 

II. POWER COUNTING SCHEMES IN NUCLEAR EFFECTIVE FIELD THE- 
ORY 

Currently, there are several power counting schemes for the nuclear effective theory 
used for multinucleon systems. Power counting means counting the power of the small 
expansion parameter of a Feynman diagram, such that one can organize the computation 
in a series expansion of this parameter. In nuclear EFT, the small expansion parameter 
Q is m^/A and p/A, where A is the cut-off scale. Here we briefly review some popular 
power counting schemes. 



In Weinberg's scheme |36l-l38| . power counting is done to the potential of the Lippmann- 
Schwinger equation, not the diagram. The leading-order (LO) potential involves the one 
pion exchange (OPE) potential and the delta function potential from contact interactions. 
Subtracting the inflnities in the LO diagram requires higher order operators with high 
power of quark mass insertions. Thus, the result has cut-off dependence that cannot 



be removed 



39 



44 1|. A similar situation happens to higher partial waves as well j41 |. 
However, within a reasonable range of cut-off, the scheme works well numerically with 
impressive flts to nucleon-nucleon (NN) scattering phase shift data at the fourth order 
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The alternative KSW scheme 



59 



60| counts the diagrams near the non-trivial UV 
flxed point of the four-nucleon operators such that the cut-off dependence is removed 
and diagrams of the same order are of equal size. The LO S-wave diagrams only contain 
non- derivative four-nucleon contact interactions, while the next-to-leading-order (NLO) 
contains OPE diagrams and diagrams with higher order four-nucleon operators. However, 



numerically, the convergence is not good in the ^5*1 channel due to the singular nature 
of the tensor pion exchan ge p otential at short distance riVj/r^, where r is the distance 
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62| . The suggests that the tensor pion exchange might not be 



between two nucleons 
perturbative. 

In view of this problem, the tensor pion exchange is resummed at the LO in the BBSvK 
scheme 40| where the ^5*0 channel follows the KSW power counting while the ^5*1 channel 
follows the Weinberg's power counting. It was shown that the cut-off can be removed in 
this scheme. 



The BKV scheme 351] seeks to fix the same problem by introducing a Pauli-Villars 
(PV) field in the ^Si channel to remove the short distance part of the singular tensor 
potential. The resulting ^Si phase shift is convergent. The price to pay is that the PV 
mass A is counted as the same order as m.^, but numerically it is close to the cut-off scale. 
However, its analytic result is very convenient to perform the matching to a pionless 
theory. Thus, we will adopt the BKV scheme in this work. 

III. THE QUARK MASS DEPENDENCE OF EFFECTIVE RANGE PARAME- 
TERS 

The S-wave nucleon-nucleon (NN) scattering amplitude is 

A=^ I , (1) 

M p cot 6 — ip^ 

where M = 938.92 MeV is the nucleon mass, p is the magnitude of the nucleon three- 
momentum in the center-of-mass (CM) frame and 6 is the S-wave phase shift. If the 



interaction (potential) is localized, then 6 has the expansion 63|, l64 1 

pcot6 = + \rop'^ + . . . , (2) 
a 2 

where the effective range parameters (ERP's) a and tq are the scattering length and the 
effective range, respectively. The shape parameter and higher order terms are not shown. 

In the BKV scheme, the amplitude can be expanded in powers of the small expansion 
parameter Q 

A = A-i + Ao + Ai + ..., (3) 
where An is of order in the expansion. Hence 



(4) 



A. Working in the BKV scheme as an explicit example 



The BKV scheme is the same as the KSW scheme in the ^5*0 channel but different in 
the ^5*1 channel The leading order (LO) amplitude of channel i arises from the diagrams 
in Fig. 5 of Ref. 
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(5) 



where Cq is the LO four-nucleon non-derivative coupling which is independent of m„. The 
next-to-leading-order (NLO) amplitude arises from the diagrams in Fig. 6 of Ref. [59j plus 
the associated diagrams with the Pauli-Villars fields 



0,6 



5^ 



0,6 



(A) 



(6) 



where e is introduced to keep track of the difference between the BKV and the KSW 
power counting, e = 1 gives the BKV result while e = gives the KSW result. 



'-'0,0 



1 + 



(/i + ip) 



+ 



1 + 
2/2 



An 



(/i + ip) 



-1 + 



4p2 



In 1 + 



4p2 



P 



ifJ' + ip) ^ 

2p 



tan 



2p 



+ 



2/2 I An 



1 



fi + ip 



+ i tan ^ 



2p 



2 



-li„ 

2 



4p2 

9 

mi 

TT 

ml + 4p2 



(7) 



where qa is the pion nucleon coupling constant, / is the pion decay constant, and we 
have imposed isospin symmetry by setting rriu = rrid and neglecting the electromagnetic 
interaction. Co,o is a NLO operator with the same structure as Co- -D2 is a non- derivative 



four-nucleon coupling with one insertion of (or m^), and C2 is a two- derivative four- 
nucleon operator that is independent of m^. /i is the renormahzation scale, and we 
have used dimensional regularization and the power-divergence subtraction procedure 
(PDS) 59| to renormalize the theory. These amplitudes are manifestly renormahzation 
scale independent order-by-order in the EFT expansion. 

Expanding the right hand side of Eq.(|l]) in powers of p, we have the matching for the 
spin singlet and triplet scattering lengths 
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where 'y^^ = ^i + An/MCf is the LO inverse scattering length. We perform the expansion 
around the physical pion mass m^'^^ ~ 138 MeV, so 7''*'' takes the physical value l/( 



phys ■ 



To fix and CqI, we just need computed at another m^^ other than m^y^. 

The matching for effective ranges gives 

rA = = ■ — h , r „„ 6 
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Unlike the scattering lengths, no lattice data is needed to study the quark mass depen- 
dence of the effective ranges since C'f' can be fixed by r^^ at m^y'^ . 

Note that the e terms in and Tq are independent, so they can be absorbed 

into counterterms Cq q and . Therefore, the KSW and BKV schemes give the same 
dependence to ERP's at NLO. 
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FIG. 1: Scattering lengths of the ^Si and ^Sq states vs. using the NLO BKV result of 
Eq.®, the physical scattering length, and the scattering length computed at m,r = 353 MeV 
with lattice QCD. The dashed(solid) lines are with(without) the higher order m.,^ dependence 
in M, / and gA included. The dot is the physical point. 



In summary, Eqs. flHllQl) can be parametrized as 
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5.423 ± 0.005 fm, rf^'^ 



0,phys 



-23. 714 ± 0.003 fm, r 
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1.764 ± 0.002 fm, a^,^, 
(0 
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2.73 ±0.03 fm) fix C2 and a combination of 7*-*^ and c/j • We 
only need a LQCD calculation of a^^^ at different to get the leading dependence 



for a^*^ and 



Currently, the smallest that a*-*) is computed on the lattice is 353.7 ±2.1 MeV 23|. 
The calculation yields a^'^^^^ = 0.63 ± 0.74 fm, a^^^"^ = 0.63 ± 0.50 fm. The central values 
yield the solid curves in Fig. 1 and 2. We can study the size of higher order corrections 
by including the dependence of M, / and qa (these are next-to-next-to-leading-order 



corrections) which is extracted from lattice data [65|, |66| to Eqs.f l8]|9|) . This yields the 
dashed curves in Fig. 1 and 2. The a^^^^^ — > 00 position can shift by ~ 20% in due to 
higher order corrections, while the corrections to a^'^o) is much smaller. When 771^^ > 100 
MeV, ^ 2 fm and is insensitive to m^. 

The analytic structure of the scattering amplitude, Eq.([l]), is that there are two cuts 
from p = zm^/2 to ioo and from p = — 2m^/2 to —ioo. There is a ^Si bound state for 
= 106 to 142 MeV (with A = 750 MeV, but the range remains the same for A = 500 to 




FIG. 2: Effective ranges of tlie and ^Sq states vs. ttItt using the NLO BKV result of Eq. 
The notations are the same as in Fig. 1. 




FIG. 3: Bound state (whenever exists) binding energy vs. for ^Si(sohd line) and ^5o (dashed 
line). The notations are the same as in Fig. 1. 

1000 MeV) and a ^5*0 bound state for = 144 to 165 MeV. The corresponding binding 
energies are shown in Fig. 3. This result can be understood by examining the scattering 
ampHtude in the effective range expansion. By keeping only the scattering length and 
effective range in Eq. ([2]) , the amphtude of Eq. ([T]) has two poles 
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If a > 0, the solution with smaller \p\ is p = ^ ^1 — yl — ^j- The bound state 
exists when < ^ < 1. On the other hand, if a < 0, the solution with smaller \p\ is 
1 - 



p 



2-ro 
a 



Since — < 0, the pole does not correspond to a bound state. 



The other pole p = ^ (l + Wl — (for both positive and negative a) is of the order of 



the ultraviolet cut-off scale 1/ro which is usually hidden in the cut starting at p = 
Thus, the bound state range is < 2ro < a, which is close to the ranges seen in Fig. 3. 
Furthermore, the maximum binding momentum is i/ro, or the maximum binding energy 
is l/(Mrg) ~ 6 MeV for tq ~ 2.5 fm. 

With the current lattice input, this theory does not have a two nucleon bound state in 
the chiral (m^ — )■ 0, or equivalently — > 0) limit. However, different power countings 
could lead to different conclusions [23]. This might indicate that m^r = 354 MeV is not 
within the common "chiral regime", i.e. within the radius of convergence of the 
expansion, for these theories. It is important to perform higher order EFT calculations 
to decide the size of the chiral regime and to answer how small m.,^ should be for future 
LQCD calculations to draw a firm conclusion about the deuteron binding energy in the 
chiral limit. 



B. A support for multiverse? 

It is curious that m^^^ is so close to the upper bound of m,r where the deuteron is 
bounded — if rn^y'^ were 5% bigger, then there would not have been deuteron at all. This 
makes it much harder for primordial nuclear synthesis to form light nuclei through the 
usual pathways and might eventually make life impossible. This interesting fine tuning 
implies that our universe sits near the edge of the parameter space where life could exist. 
In Ref. [g^], it is argued that this is not a fine tuning but a natural case if multiverse exists: 
In a multiverse, which is an ensemble of many universes including ours, the majority of the 
universes do not allow life to exist since it requires lots of conditions to be satisfied. Thus, 
the peak of the distribution in this multiverse will be more likely to sit outside the 
parameter space where life is possible. In that case, the tail of the distribution goes across 
this parameter space and then one finds that most of the universes that permits life is 
near the edge of the parameter space. Thus, if the multiverse exists, without fine tuning, 

le edge of the parameter space where hfe is possible (called 



our universe should live near t 

the catastrophic boundary in [67|). It is interesting to note that our case of the deuteron 
bound state is consistent with this pattern, similar to the example of the cosmological 
constant whose value is close to the allowed range obtained by Weinberg through the 



anthropic principle 68|, l69|] and several examples worked out in 67| . Although we do not 



consider this as a sharp test of the multiverse conjecture, because the conjecture cannot be 



falsified even if the physical is far away from the edge of the allowed parameter space, 
it is still interesting to see whether there are cases being consistent with this conjecture. 



IV. MATCHING BETWEEN THE THEORY WITH AND WITHOUT PIONS 

We are interested in using the theory without pions to describe low energy processes 
(where p < m.^ so pions can be integrated out) at non-physical m.^. The matching between 
the pionful and pionless EFT's at those gives the m^r dependence of the couplings in 
the pionless theory. Those couplings are the ERP's of NN scattering mentioned above 
and current operators when coupled to external currents. 

We can classify the non-derivative single-nucleon (one-body) current operators by how 
they transform in the spin-isospin space: the scalar-scalar operator (N'^N), scalar- vector 
operator (A^VjA^), vector-scalar operator (A^VjA^) and vector- vector operator (N'^aiTjN), 
where o"j(rj) acts on the spin(isospin) space and the spacial indexes i,j = 1,2,3. The 
non-derivative scalar-scalar and scalar- vector operators originate from matrix elements of 
the quark level operators g7og and g7oTjg. They do not have two-body currents due to 
vector current conservation. For vector-scalar currents, they could originate from matrix 
elements of the isoscalar quark axial operator qjij^q or the magnetic part of the vector 
current q'^iQ, so the corresponding two body-currents exist. For vector-vector currents, the 
corresponding quark level operator is q'Ji'J^Tjq and the two body-currents (called Gamow- 
Teller operators) also exist. 

From matching the isoscalar magnetic current between the theory with and without 
pions, we conclude that the vector-scalar two-body currents do not depend on pion mass 



at the leading order [70|]. The matching of the two-body Gamow- Teller operator |7ll. |72| 
yields 



2772 °V^^ + 27 

ea^'^^o^m^ + ml {^a^^^^'^m^ - 4 j 7 - 2771^7^ (^a^'^^^m^ - 5 



Qa'^pm'l (m^ + 27) 

-27^ fsa^''^''^^ + e) + 12a('^o)7^1 , (12) 



where Iqt is m.,, independent, 7 = ^1 — 1 — 2rQ^^^ / a'^'^^^^ /r^ '^^^ is the deuteron bind- 
ing momentum and Hi is the single nucleon coupling (for the isovector magnetic current, 
Ki is the isovector nucleon magnetic moment; for weak coupling, ki is proportional to 
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FIG. 4: Deuteron charge radius {\J ij""^)-, magnetic moment (fJ-Ai), and electric polarizability 
(oeo) vs. m^. The dots are the physical points. 

qa)- There is no unknown parameter in the m,r dependent term. 

V. THE QUARK MASS DEPENDENCE OF MORE TWO NUCLEON OBSERV- 
ABLES 

In this section we apply the m.^ dependent couphngs in the pionless EFT, which has 
been worked out in the previous sections, to compute several physical observables involving 
deuterons. 



A. Deuteron properties 

The deuteron charge radius has the expression 



33] 



72 (1 - 7prf) 



(13) 



where the isoscalar charge radius of the nucleon {i^n,o) = 0.79 ±0.01 fm and pd 
As expected, the deuteron charge radius (r^) is set by the inverse binding momentum 
1/7 when the nucleon charge radius is neghgible. The dependence for = 125 — 141 



MeV (where deuteron is bounded) is s^ 
The deuteron magnetic moment is 



lown in Fig. 4. 
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where = 0.44 is the nucleon isoscalar magnetic moment in units of nuclear magneton 
(N.M.) and the vector-scalar two-body current Ly^s is "^tt independent. Neither the 
one- nucleon nor the two-nucleon contribution is sensitive to m^. The sum is also shown 
in Fig. 4. 



The deuteron polarizability is computed as 



33| 



aM 



(15) 



3274 (1-7P,)' 

a = 1/137 is the fine structure constant. It has a strong m,r dependence as is shown in 
Fig. 4. 



B. Reaction: np o d-y 



The process np -H- d'y is relevant for BBN. Its cross section is proportional to the wave 
function overlap between the initial and final states. Since the deuteron size I/7 is very 
sensitive to near rn^y^ , the cross section also changes dramatically in this region. 

The total cross section for np — )■ (^7 is 



a {np — )■ d-y) 



Ana (72 + 



-f'^M^p 

where p is the magnitude of the momentum of each nucleon in the center-of-mass frame. 
The electric dipole {El) transition yields 

p^W^-i^ 



\x 



El\ 



(72 + p"^) 

The Magnetic dipole (Ml) transition yields 
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FIG. 5: The cross section for jd — ?> np as a function of the incident photon energy in MeV. The 
sohd curve is for the physical m^r, while the dotted, dotdashed, and dashed curves are for 771,7^ = 
125, 130, and 141 MeV, respectively. 

The Gamow- Teller two-body current: Lgt = —4.513 fm^ at = ■m^'^y^, is fitted from 



the measured cross section a'^^P* = 334.2 ± 0.5 mb 73j] using incident neutrons of speed 
\v\ = 2200 m/s. The dependence of Lgt is shown in Eq.([2]). The isovector nucleon 
magnetic moment /ti = 2.35 — 5'^M(m^ — m^^y^)/ (27r/^), where we have applied the m.,^ 



dependence calculated from ChPT 2l(]. The cross section of the reverse process with 
deuteron being at rest is 

a {-id np) = a {np d-f) , (17) 

where is the incident photon energy. This deuteron photo-disintegration cross section 
for = 125 — 141 MeV is shown in Fig. 5. 



VI. CONCLUSION 



We have studied the implications of lattice QCD determinations of the S-wave nucleon- 
nucleon scattering lengths at unphysical light quark masses. It is found that with the help 
of nuclear effective field theory, not only the quark mass dependence of the effective range 
parameters, but also the leading quark mass dependence of all the low energy (with 
p ^ rriTr) deuteron matrix elements can be obtained. The quark mass dependence of 



deuteron charge radius, magnetic moment, polarizability and the deuteron photodisinte- 
gration cross section are shown based on the NPLQCD lattice calculation of the scattering 
lengths at 354 MeV pion mass and the NEFT power counting scheme of Beane, Kaplan 
and Vuorinen. Further improvement can be obtained by performing the lattice calculation 
at smaller quark masses. But at the same time, it is important to perform higher order 
EFT calculations to decide the radius of convergence in in order to answer how small 
should be for future LQCD calculations to provide reliable m.^ dependence for two 
nucleon observables all the way to the chiral limit. 

Our result can be used to constrain the time variation of isoscalar combination of u 
and d quark mass ruq, to help the anthropic principle study to find the nig range which 
allows the existence of life, and to provide a weak test of the multiverse conjecture. 

We thank Martin Savage for providing us the formulae describing the LQCD m„ de- 
pendence of M, / and qa- This work is supported by the NSC and NCTS of ROC. 
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